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1. Introduction 

Group relaxations of integer programs were introduced by Ralph Gomory in the 1960s 
|13]| , [|15|] . Given a general integer program of the form 

(1) minimize {c ■ x : Ax = b, x > 0, integer}, 

its group relaxation is obtained by dropping non-negativity restrictions on all the basic 
variables in the optimal solution of its linear relaxation. In this paper, we survey recent 
results on group relaxations obtained from the algebraic study of integer programming using 
Grobner bases of toric ideals |[36|| . No knowledge of these methods is assumed, and the 
exposition is self-contained and hopefully accessible to a person familiar with the traditional 
methods of integer programming. Periodic comments on the algebraic origins, motivations 
and counterparts of many of the described results — which the reader may pursue if desired 
— offer a more complete picture of the theory. 

For the sake of brevity, we will bypass a detailed account of the classical theory of group 
relaxations. A short expository account can be found in |33], §24.2], and a detailed set of 
lecture notes on this topic in p^] . We give a brief synopsis of the essentials based on the 
recent survey article by Aardal et. al and refer the reader to any of the above sources for 
further details and references on the classical theory of group relaxations. 

Assuming that all data in ([!]) are integral and that As is the optimal basis of the linear 
relaxation of ([!]), Gomory's group relaxation of (Q) is the problem 



(2) 



minimize {c ■ xn '■ A^A^xn = A^b {mod 1), xn > 0, integer}. 
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Here B and iV are the index sets for the basic and non-basic columns of A corresponding 
to the optimal solution of the linear relaxation of (HI). The vector xn denotes the non-basic 
variables and the cost vector c = — cbA]^A n where c = (c#, c^) is partitioned according 
to B and N. The notation A^A^xn = A^b (mod 1) indicates that A^AnXn — A B b is a 
vector of integers. Problem @) is called a "group relaxation" of ([!]) since it can be written 
in the canonical form 

(3) minimize {c ■ Xn : 'S^gjXj = go (mod G), x^ > 0, integer} 

j&N 

where G is a finite abelian group and gj G G. Problem (0) can be viewed as a shortest 
path problem in a graph on \G\ nodes which immediately furnishes algorithms for solving 
it. Once the optimal solution x* N of (0) is found, it can be uniquely lifted to a vector 
x* = (x* B ,x* N ) G Z7 1 such that Ax* = b. If x* B > then x* is the optimal solution of (|l|). 
Otherwise, c- x* is a lower bound for the optimal value of (|l|). Several strategies are possible 
when the group relaxation fails to solve the integer program. See ||,[I7|, [31 and jKJ for 



work in this direction. A particular idea due to Wolsey ||2] that is very relevant for this 
paper is to consider the extended group relaxations of ([!]). These are all the possible group 
relaxations of ([I]) obtained by dropping non-negativity restrictions on all possible subsets 
of the basic variables xb in the optimum of the linear relaxation of (|I|). Gomory's group 
relaxation (§) of ([!]) and ([[]) itself are therefore among these extended group relaxations. 
If (H) does not solve (|l|), then one could resort to other extended relaxations to solve the 
problem. At least one of these extended group relaxations (in the worst case (0) itself) is 
guaranteed to solve the integer program ([!]). 



The convex hull of the feasible solutions to (g) is called the corner polyhedron [14]. A 
major focus of Gomory and others who worked on group relaxations was to understand the 
polyhedral structure of the corner polyhedron. This was achieved via the master polyhedron 



of the group C [15] which is the convex hull of the set of points 



{z : ^^^gzg = g (mod G), z > 0, integer}. 
g eG 

Facet-defining inequalities for the master polyhedron provide facet inequalities of the corner 
polyhedron |T3. As remarked in H], this landmark paper |T3[ introduced several of the now 



standard ideas in polyhedral combinatorics like projection onto faces, subadditivity, master 
polytopes, using automorphisms to generate one facet from another, lifting techniques and 
so on. See |TE[ for further results on generating facet inequalities. 



In the algebraic approach to integer programming, one considers the entire family of 
integer programs of the form (0) as the right hand side vector b varies. Definition |2.6| defines 
a set of group relaxations for each program in this family. Each relaxation is indexed by 
a face of a simplicial complex called a regular triangulation (Definition [2.1|) . This complex 



encodes all the optimal bases of the linear programs arising from the coefficient matrix A 
and cost vector c (Lemma |2.3|) . The main result of Section 2 is Theorem |2.8| which states 



that the group relaxations in Definition |2.6| are precisely all the bounded group relaxations 
of all programs in the family. In particular, they include all the extended group relaxations 
of all programs in the family and typically contain more relaxations for each program. This 
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theorem is proved via a particular reformulation of group relaxations which is crucial for the 
rest of the paper. This and other reformulations are described in Section 2. 

The most useful group relaxations of an integer program are the "least strict" ones among 
all those that solve the program. By this we mean that any further relaxation of non- 
negativity restrictions will result in group relaxations that do not solve the problem. The 
faces of the regular triangulation indexing all these special relaxations for all programs in the 
family are called the associated sets of the family (Definition |Q| ). In Section 3 we develop 
tools to study associated sets. This leads to Theorem |3.10| which characterizes associated 



sets in terms of standard pairs and standard polytopes. Theorem |3.11| shows that one can 
"read off" the "least strict" group relaxations that solve a given integer program in the 
family from these standard pairs. 

The results in Section 3 lead to an important invariant of the family of integer programs 
being studied called its arithmetic degree. In Section 4 we discuss the relevance of this 



invariant and give a bound for it based on a result of Ravi Kannan (Theorem 4.8). His 
result builds a bridge between our methods and those of Kannan, Lenstra, Lovasz, Scarf and 
others that use geometry of numbers in integer programming. 

Section 5 examines the structure of the poset of associated sets. The main result in 



this section is the chain theorem (Theorem 5.2 ) which shows that associated sets occur in 



saturated chains. Theorem ^]4| bounds the length of a maximal chain. 

In Section 6 we define a particular family of integer programs called a Gomory family, 
for which all associated sets are maximal faces of the regular triangulation. Theorem ^]2| 
gives several characterizations of Gomory families. We show that this notion generalizes the 



classical notion of total dual integrality in integer programming []33], §22]. We conclude in 
Section 7 with constructions of Gomory families from matrices whose columns form a Hilbert 
basis. In particular, we recast the existence of a Gomory family as a Hilbert cover problem. 



This builds a connection to the work of Sebo |[34|| , Bruns & Gubeladze J7J and Firla & Ziegler 



TT| on Hilbert partitions and covers of polyhedral cones. We describe the notions of super 



and A-normality both of which give rise to Gomory families (Theorems [?\7] and |7.14p . 



The majority of the material in this paper is a translation of algebraic results from [^T 



22 1, P3J, PR §8 and §12. D], |BB| and [35]. The translation has sometimes required new 



definitions and proofs. Kannan's theorem in Section 4 has not appeared elsewhere. 

We will use the letter N to denote the set of non-negative integers, WL to denote the real 
numbers and Z for the integers. The symbol P C Q denotes that P is a subset of Q, possibly 
equal to Q, while P C Q denotes that P is a proper subset of Q. 

2. Group Relaxations 

Throughout this paper, we fix a matrix A G Z dxn of rank d, a cost vector c G Z n and 
consider the family IPa,c of all integer programs 

IPa,c{o) '■= minimize {c ■ x : Ax = b, x G N n } 

as b varies in the semigroup NA := {Au : u G N n } C Z d . This family is precisely the set of 
all feasible integer programs with coefficient matrix A and cost vector c. The semigroup NA 
lies in the intersection of the <i-dimensional polyhedral cone cone(A) := {Au : u > 0} C M. d 
and the <i-dimensional lattice 7*A := {Au : u G Z™} C Z d . For simplicity, we will assume 
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that cone(A) is pointed and that {mgP : Au = 0}, the kernel of A, intersects the non- 
negative orthant of M. n only at the origin. This guarantees that all programs in IPa, c are 
bounded. In addition, the cost vector c will be assumed to be generic in the sense that each 
program in IPa, c has a unique optimal solution. 

The linear relaxation of IPa, c {^) is the linear program 

LPa, c (o) '■= minimize { c ■ x : Ax = b, x > 0}. 

We denote by LPa, c the family of all linear programs of the form LPa, c {o) as b varies in 
cone(A). These are all the feasible linear programs with coefficient matrix A and cost vector 
c. Since all data are integral and all programs in IPa, c are bounded, all programs in LPa, c 
are bounded as well. 

In the classical definitions of group relaxations of IPA. c {b), one assumes knowledge of the 
optimal basis of the linear relaxation LPA. c {b). In the algebraic set up, we define group 
relaxations for all members of IPa, c at one shot and, analogously to the classical setting, 
assume that the optimal bases of all programs in LPa jC are known. This information is 
carried by a polyhedral complex called the regular triangulation of cone(A) with respect to c. 

A polyhedral complex A is a collection of polyhedra called cells (or faces) of A such that: 

(i) every face of a cell of A is again a cell of A and, 

(ii) the intersection of any two cells of A is a common face of both. 

The set-theoretic union of the cells of A is called the support of A. If A is not empty, then 
the empty set is a cell of A since it is a face of every polyhedron. If all the faces of A are 
cones, we call A a cone complex. 

For a C {1, . . . , n}, let A a be the sub matrix of A whose columns are indexed by a, and 
let cone(A a ) denote the cone generated by the columns of A a . The regular subdivision A c of 
cone(A) is a cone complex with support cone(A) defined as follows. 

Definition 2.1. For a C {1,... ,n}, cone(A a ) is a face of the regular subdivision A c of 
cone(A) if and only if there exists a vector y G M. d such that y ■ aj = Cj for all j 6 a and 
y ■ dj < Cj for all j ^ a. 

The regular subdivision A c can be constructed geometrically as follows. Consider the cone 
in M. d+1 generated by the lifted vectors G IR d+1 where a« is the ith column of A and 

Ci is the ith component of c. The lower facets of this lifted cone are all those facets whose 
normal vectors have a negative (d + l)th component. Projecting these lower facets back onto 
cone(A) induces the regular subdivision A c of cone(A) (see JJ). Note that if the columns 
of A span an affine hyperplane in M d , then A c can also be seen as a subdivision of conv(A), 
the (d — l)-dimensional convex hull of the columns of A. 

The genericity assumption on c implies that A c is in fact a triangulation of cone(A) (see 
37j). We call A c the regular triangulation of cone(A) with respect to c. For brevity, we 



may also refer to A c as the regular triangulation of A with respect to c. Using a to label 
cone(A a ), A c is usually denoted as a set of subsets of {1, . . . , n}. Since A c is a complex of 
simplicial cones, it suffices to list just the maximal elements (with respect to inclusion) in 
this set of sets. By definition, every one dimensional face of A c is of the form cone(aj) for 
some column aj of A. However, not all cones of the form cone(ai), a« a column of A, need 
appear as a one dimensional cell of A c . 



AN ALGEBRAIC PERSPECTIVE OF GROUP RELAXATIONS 

w s'is' 




4 • 
3 O 
2 O 



(c) 



Figure 1. Regular triangulation A c for c = (1, 0, 0, 1) (Example O (i)). 
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Figure 2. Regular triangulation A c / for c' = (0, 1, 0, 1) (Example |2T^ (ii)). 



Example 2.2. (i) Let A 



and c = (1, 0, 0, 1). The four columns of A are 



4. Figure p] (a) shows 



1111 
12 3 

the four dark points in Figure [I] labeled by their column indices 1 

the cone generated by the lifted vectors (a*, Cj) e M 3 . The rays generated by the lifted vectors 
have the same labels as the points that were lifted. Projecting the lower facets of this lifted 
cone back onto cone(A), we get the regular triangulation A c of cone(A) shown in Figure [I] 
(b). The same triangulation is shown as a triangulation of conv(A) in Figure [I] (c). The 
faces of the triangulation A c are {1, 2}, {2, 3}, {3, 4}, {1}, {2}, {3}, {4} and 0. Using only the 
maximal faces, we may write A c = {{1, 2}, {2, 3}, {3, 4}}. 

(ii) For the A in (i), cone(A) has four distinct regular triangulations as c varies. For instance, 
the cost vector d = (0, 1, 0, 1) induces the regular triangulation A c / = {{1, 3}, {3, 4}} shown 
in Figure @ (b) and (c). Notice that {2} is not a face of A c /. 

1 ' 2 1 and c= (1,0, 0,1), then A c = {{1,2}, {2, 3}, {3, 4}}. However, in 



(m)l£A- x 2 3 
this case, A c can only be seen as a triangulation of cone(A) and not of conv(A) 



□ 
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For a vector x G M. n , let supp(x) := {i : Xi ^ 0} denote the support of x. The significance 
of regular triangulations for linear programming is summarized in the following proposition. 

Proposition 2.3. [[37], Lemma 1.4] An optimal solution of LPa, c {o) is any feasible solution 
x* such that supp(x*) = r where r is the smallest face of the regular triangulation A c such 
that b G cone(A T ). 

Proposition |2l| implies that o C {1, . . . , n} is a maximal face of A c if and only if A a 
is an optimal basis for all LPA }C (b) with b in cone(A a ). For instance, in Example [2.2| (i), 
if b = (4, 1)* then the optimal basis of LPA )C {b) is [01,02] where as if b = (2,2)*, then the 
optimal solution of LPa iC (6) is degenerate and either [ai,a 2 ] or [a 2 , 03] could be the optimal 
basis of the linear program. (Recall that a; is the ith column of A.) All programs in LPa, c 
have one of [ai, a 2 ], [a 2 , a 3 ] or [a 3 , a 4 ] as its optimal basis. 

Given a polyhedron P C W 1 and a face F of P, the normal cone of F at P is the cone 
Np(F) := {u £ W 1 : a; • a;' > uj ■ x, ior all x' & F and x G P}. The normal cones of all faces 
of P form a cone complex in M n called the normal fan of P. 

Proposition 2.4. TTie regular triangulation A c of cone(A) is the normal fan of the polyhe- 
dron P c := {y G W 1 :yA< c}. 

Proof. The polyhedron P c is the feasible region of maximize {y ■ b : yA < c, y G the 
dual program to LP^ >C (6). The support of the normal fan of P c is cone(A), since this is the 
polar cone of the recession cone {y e R d : yA < 0} of P c . Suppose b is any vector in the 
interior of a maximal face cone(A a ) of A c . Then by Proposition |2.3| , LPA tC {b) has an optimal 
solution x* with support a. By complementary slackness, the optimal solution y to the dual 
of LPA, c {b) satisfies y ■ aj = Cj for all j G a and y ■ aj < Cj otherwise. Since a is a maximal 
face of A c , y-dj < cj for all j SjL a. Thus y is unique, and cone{A a ) is contained in the normal 
cone of P c at the vertex y. If b lies in the interior of another maximal face cone(A T ) then y ; , 
(the dual optimal solution to LPa )C (6)) satisfies y' ■ A T = c T and y' ■ A f < c f where r 7^ a. As 
a result, y' is distinct from 7/, and each maximal cone in A c lies in a distinct maximal cone 
in the normal fan of P c . Since A c and the normal fan of P c are both cone complexes with 
the same support, they must therefore coincide. □ 



Example |2.2| continued. 



Figure |3] (a) shows the polyhedron P c for Example |2.2j (i) with all its normal cones. The 
normal fan of P c is drawn in Figure |3| (b). Compare this fan with that in Figure [I] (b). □ 



Corollary 2.5. The polyhedron P c is simple if and only if the regular subdivision A c is a 
triangulation ofcone(A). 

Regular triangulations were introduced by Gel'fand, Kapranov and Zelevinsky and 
have various applications. They have played a central role in the algebraic study of integer 
programming (|[36||, |37|), and we use them now to define group relaxations of IPa, c (P)- 

A subset r of {1, . . . , n} partitions and x f where x T consists of the 

variables indexed by r and x f the variables indexed by the complementary set f. Similarly, 
the matrix A is partitioned as A = [A T , A f ] and the cost vector ). If a is a 

maximal face of A c , then A a is nonsingular and Ax = b can be written as x a = y4 o : 1 (6— A^x^). 
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Figure 3. The polyhedron P c and its normal fan for Example |2.2| (i). 

Then c-x = c a {A~ x {b- A^Xe )) + c^x a = c a A~ x b J r {c^-c a A~ x A^)x^. Let c 5 := c^-c^A^A^ 
and, for any face r of a, let Cf be the extension of c s to a vector in M' T ' by adding zeros. 
We now define a group relaxation of IPa, c (P) with respect to each face r of A c . 

Definition 2.6. The group relaxation of the integer program IPa., c (P) with respect to the 
face t of A c is the program: 

G T (b) = minimize {c f ■ x f : A T x T + A f x f = b, x f > 0, (x T , x f ) G Z n }. 

Equivalently, G T (b) = minimize {c f -x f : A f x f = b (mod 7*A T ), x f > 0, integer} where 
r LA r is the lattice generated by the columns of A T . Suppose x% is an optimal solution to the 
latter formulation. Since r is a face of A c , the columns of A T are linearly independent, and 
therefore the linear system A T x T + A f x% = b has a unique solution. Solving this system for 
x T , the optimal solution x% of G T {b) can be uniquely lifted to the solution (x* , x* f ) of Ax = b. 
The formulation of G T (b) in Definition [2.6| shows that x* T is an integer vector. The group 
relaxation G T (b) solves /Pa iC (6) if and only if x* is also non- negative. 

The group relaxations of IPA ;C {b) from Definition ETBI contain among them the classical 
group relaxations of IPA, c {b) found in the literature. The program G a (b), where A a is the 
optimal basis of the linear relaxation LPA tC {b), is precisely Gomory's group relaxation of 
IPA,c{b) HI]. The set of relaxations G T (b) as r varies among the subsets of this a are the 



extended group relaxations of IPA, c {b) defined by Wolsey pL2| . Since G A c , G®(b) = IPA :C (b) 



is a group relaxation of /Pa c (6), and hence /Pa c (6) will certainly be solved by one of its 



8 



REKHA R. THOMAS 



extended group relaxations. However, it is possible to construct examples where a group 
relaxation G T (b) solves JPa jC (6), but G T (b) is neither Gomory's group relaxation of IPa, c (P) 
nor one of its nontrivial extended Wolsey relaxations (see Example ^4.2| ). Thus, Definition [2.6| 
typically creates more group relaxations for each program in IPa, c than in the classical 
situation. This has the obvious advantage that it increases the chance that JP^ )C (6) will 
be solved by some non-trivial relaxation, although one may have to keep track of many 



more relaxations for each program. In Theorem |2.8| , we will prove that Definition |2.6| is the 
best possible in the sense that the relaxations of IPA, c {b) defined there are precisely all the 
bounded group relaxations of the program. 

The goal in the rest of this section is to describe a useful reformulation of the group 
problem G T (b) which is needed in the rest of the paper and in the proof of Theorem p7S 



Given a sublattice A of Z™, a cost vector w E W 1 and a vector v E N n , the lattice program 
defined by this data is 

minimize {w ■ x : x = v (mod A), x E N n }. 

Let C denote the (n — <i)-dimensional saturated lattice {x E Z n : Ax = 0} C Z n and u 
be a feasible solution of the integer program IPA, c (b)- Since IPA, c (b) = minimize {c ■ x : 
Ax = b (= Au), x E N 11 } can be rewritten as minimize {c ■ x : x — u E C, x E N n }, 
IPa,c(P) is equivalent to the lattice program 

minimize {c ■ x : x = u (mod £), x E M n }. 

For r C {1, . . . , n}, let ir T be the projection map from R" — > M' r ' that kills all coordinates 
indexed by r. Then C T := 7r r (£) is a sublattice of Zl T l that is isomorphic to C: Clearly, n T : 
£ — > C T is a surjection. If tt t (v) = tt t (v') for v, v' E C, then A T v T + A f v f = = A T v' T + A f v' f , 
implies that A T (v T — v' T ) = 0. Then v T = v' T since the columns of A T are linearly independent. 
Using this fact, G T (b) can also be reformulated as a lattice program: 



G T (b) = minimize {c f ■ x f 

= minimize {c f ■ x f 

= minimize {c f ■ x f 

= minimize {c f ■ x f 



A t x t ~\~ AfX^- — 6, Xf ^ 0, (x r , Xf^j E Z^j 
(x T ,x f y - (u T ,u f y E £, x f E N |f| } 
Xf — Uf E £ T , Xf E N |f| } 
Xf = n T (u) (mod C T ), Xf E N' r '}. 

Lattice programs were shown to be solved by Grobner bases in []39|j . Theorem 5.3 in 
gives a geometric interpretation of these Grobner bases in terms of corner polyhedra. 
This paper was the first to make a connection between the theory of group relaxations and 
commutative algebra (see [^, §6]). Special results are possible when the sublattice A is of 



finite index. In particular, the associated Grobner bases are easier to compute. 

Since the (n — d)-dimensional lattice C C Z n is isomorphic to C a C Z' "' for o~ E A c , C a is 
of finite index if and only if a is a maximal face of A c . Hence the group relaxations G a (b) 
as a varies over the maximal faces of A c are the easiest to solve among all group relaxations 
of IPa,c{V)- They contain among them Gomory's group relaxation of IPa, c - We give them a 
collective name in the following definition. 

Definition 2.7. The group relaxations G a (b) of IPA, c (b), as a varies among the maximal 
faces of A c , are called the Gomory relaxations of IPa, c (P)- 
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It is useful to reformulate G T (b) once again as follows. Let B G Z nx( - n ~ d ^ be any matrix 
such that the columns of B generate the lattice £, and let u be a feasible solution of IPa, c (P) 
as before. Then 

IPa,c(P) — minimize {c ■ x : x — u G C, x E N n } 

= minimize {c ■ x : x = u — Bz, x > 0, z6 r L n ~ d }. 
The last problem is equivalent to minimize {c ■ (u — Bz) : Bz < u, z G Z n ~ d } and, 
therefore IPA, c {b) is equivalent to the problem 

(4) minimize {{-cB) ■ z : Bz < u, z E Z n ~ d }. 

There is a bijection between the set of feasible solutions of (^) and the set of feasible solutions 
of IPj^ c (b) via the isomorphism z i— > u — Bz. In particular, G M. n ~ d is feasible for (f|) and 
it is the pre-image of u under this map. 

If B T denotes the |f| x (n — d) submatrix of B obtained by deleting the rows indexed by 
r, then C T = n T (C) = {B T z : z G Z n ~ d }. Using the same techniques as above, G T (b) can 
be reformulated as 

minimize {(—c f B T ) ■ z : B T z < n T (u), z G Z n ~ d }. 

Since c f = 7r T (c — c a A~ l A) for any maximal face a of A c containing r and the support of 
c — c a A~ x A is contained in f , Cf-B T = (c — c a A~ 1 A)B = cB since = 0. Hence G T (b) is 
equivalent to 

(5) minimize {{-cB) ■ z : B f z < ir T (u), z G Z n ~ d }. 

The feasible solutions to are the lattice points in the rational polyhedron P u := {z G 
IR n ~ d : < u}, and the feasible solutions to @ are the lattice points in the relaxation 
:= {z G IR n_d : B T z < n T (u)} of P u obtained by deleting the inequalities indexed by r. 
In theory, one could define group relaxations of IPA, c {b) with respect to any r C {1, . . . , n}. 
The following theorem illustrates the completeness of Definition |2.6| . 

Theorem 2.8. The group relaxation G T (b) of JPa )C (6) /ias a finite optimal solution if and 
only if t Q {1, . . . ,n} is a face of A c . 

Proof. Since all data are integral it suffices to prove that the linear relaxation 

minimize {(—cB) ■ z : z G P^} 

is bounded if and only if r G A c . 

If t is a face of A c then there exists y G M. d such that yA T = c T and yA f < c f . Using the fact 
that A T B T + A f B f = we see that cB = c T B T +c f B f = yA T B T + c f B f = y(-A f B f ) + c f B f = 
(c f — yA f )B T . This implies that cB is a positive linear combination of the rows of B T since 
Cf — yA f > 0. Hence cB lies in the polar of {z G IR n ~ d : B T z < 0} which is the recession 
cone of PJ" proving that the linear program minimize {(— cB) ■ z : z G P^} is bounded. 

The linear program minimize {(— cP) ■ z : 2 G P^} is feasible since is a feasible 
solution. If it is bounded as well then minimize {c T x T + c f x f : A T x T + A f x f = b, x f > 0} 
is feasible and bounded. As a result, the dual of the latter program maximize {y-b : yA T = 
c T , yAf < Cf} is feasible. This shows that a superset of r is a face of A c which implies that 
t G A c since A c is a triangulation. □ 
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3. Associated Sets 

The group relaxation G T (b) (seen as (|D) solves the integer program IPa, c (P) (seen as (f|)) 
if and only if both programs have the same optimal solution z* £ Z n ~ d . If G T {b) solves 
IPA, c {b) then G T (b) also solves JPa iC (6) for every r' C r since G T (b) is a stricter relaxation 
of IPA )C (b) than G T {b). For the same reason, one would expect that G T {b) is easier to solve 
than G T (b). Therefore, the most useful group relaxations of IPA, c (b) are those indexed by 
the maximal elements in the subcomplex of A c consisting of all faces r such that G T (b) solves 
IPA, c {b)- The following definition isolates such relaxations. 

Definition 3.1. A face r of the regular triangulation A c is an associated set of IPa, c (or 

is associated to IPa,c) if f° r some b £ NA, G T {b) solves IP A , c ( b ) but G ( b ) does not for a11 
faces t' of A c such that r C r'. 

The associated sets of IPa, c carry all the information about all the group relaxations 
needed to solve the programs in IPa, c - in this section we will develop tools to understand 
these sets. We start by considering the set O c C N n of all the optimal solutions of all 
programs in IPa, c - A basic result in the algebraic study of integer programming is that O c 
is an order ideal or down set in N n , i.e., if u £ O c and v < u, v G N n , then v G O c . One way 
to prove this is to show that the complement M c := N n \(9 c has the property that if v G M c 
then v + W 1 C M c . Every lattice point in W 1 is a feasible solution to a unique program in 
IPa,c (u £ N n is feasible for IPj^ c {Auf). Hence, M c is the set of all non-optimal solutions of 
all programs in IPa, c - A set P C N n with the property that p + N"CP whenever p G P 
has a finite set of minimal elements. Hence there exists a%, . . . , at G Af c such that 

t 

K = {J(a t + N n ). 
i=i 

As a result, O c is completely specified by the finitely many "generators" a±, ... ,a t of its 
complement M c . See [^0j for proofs of these assertions. 

Example 3.2. Consider the family of knapsack problems: 

minimize {lOOOCtei + 100a;2 + £3 : 2xi + 5x 2 + 8x 3 = b, (xi,x 2 , x 3 ) G N 3 } 

as b varies in the semigroup N[2 5 8]. The set Af c is generated by the vectors 

(0, 8, 0), (1, 0, 1), (1, 6, 0), (2, 4, 0), (3, 2, 0), and (4, 0, 0) 

which means that Af c = ((0, 8, 0) + N 3 ) U • • ■ U ((4, 0, 0) + N 3 ). Figure g is a picture of Af c 
(created by Ezra Miller). The white points are its generators. One can see that O c consists 
of finitely many points of the form (p, q, 0) where p > 1 and the eight "lattice lines" of points 
(0,i,*),i = 0,... ,7. □ 

For the purpose of computations, it is most effective to think of J\f c and O c algebraically. 
A monomial x u in the polynomial ring S := k[x±, . . . ,x n ] is a product x u = x^x^ 2 ■ ■ -x^™ 
where u = (ui, . . . , u n ) G N n . We assume that k is a field, say the set of rational numbers. 
For a scalar k u G k and a monomial x u in S, we call k u x u a term of S. A polynomial 
f = k u x u in S is a combination of finitely many terms in S. A subset / of S is an ideal 
of S if (1) I is closed under addition, i.e., f,g £ I f + g £ I and (2) if f £ I and 
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g G S then fg G I. We say that / is generated by the polynomials fi, ■ ■ ■ , ft, denoted as 
I = (fx, . . . , f t ), if I = {Y^i=i fi9i '■ 9i By Hilbert's basis theorem, every ideal in S 

has a finite generating set. An ideal M in S is called a monomial ideal if it is generated by 
monomials, i.e., M = (x Vl , ... ,x vt ) for monomials x Vl , ... , x Vt in S 1 . The monomials that do 
not lie in M are called the standard monomials of M. The cosi of a term /c u x u with respect 
to a vector c G R n is the dot product c-u. The initial term of a polynomial f — Yl k u x u G 5 
with respect to c, denoted as in c (f), is the sum of all terms in / of maximal cost. For any 
ideal I C S, the initial ideal of I with respect to c, denoted as in c (I), is the ideal generated 
by all the initial terms in c (f) of all polynomials / in I. These concepts come from the theory 
of Grobner bases for polynomial ideals. See for an introduction. 

The toric ideal of the matrix A, denoted as I a, is the binomial ideal in S defined as: 

I A ■= ( x u - x v : M,ueN" and Au = Av). 

Toric ideals provide the link between integer programming and Grobner basis theory. See 
|36| and [fy] for an introduction to this area of research. This connection yields the following 
basic facts that we state without proofs. (Recall that the cost vector c of IPa, c was assumed 
to be generic in the sense that each program in IPa.c has a unique optimal solution.) 

Lemma 3.3. (i) If c is generic, then the initial ideal in c (I a) is a monomial ideal, 
(ii) A lattice point u is non-optimal for the integer program IPa, c (Au), or equivalently, 
u G Af c , if and only if x u lies in the initial ideal in c (lA)- In other words, a lattice point u 
lies in O c if and only if x u is a standard monomial of in c (lA)- 

(Hi) The reduced Grobner basis Q c of I a with respect to c is the unique minimal test set 
for the family of integer programs IPa,c- 
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(iv) If u is a feasible solution of IPa, c (P); and x u * is the unique normal form of x u with 
respect to Q c , then u* is the optimal solution of IPa, c {o)- 

We do not elaborate on parts (iii) and (iv) of Lemma |3.3| . They are not needed for what 
follows and are included for completeness. Since c is generic, Lemma |3.3| (ii) implies that 
there is a bijection between the lattice points of O c and the semigroup NA via the map 
(ft a '■ O c — > NA such that u \— > Au. The inverse of (ft a sends a vector b G NA to the optimal 
solution of IPA, c {b). 



Example |3.2| continued. In this example, the toric ideal I a = (x\ — £3, x\ — x\x%) and its 



initial ideal with respect to the cost vector c = (10000, 100, 1) is 

^c(/a) = (^2' X\Xj,, X1X2, X1X2, X1X2, X±). 

Note that the exponent vectors of the generators of in c {lA) are the generators of M c - d 

We will now describe a certain decomposition of the set O c which in turn will shed light on 
the associated sets of IPa, c - For u £ N n , consider Q u := {z G M n_d : Bz < u, (-cB)-z < 0} 
and its relaxation Q f u := {z G M n " d : B f z < tt t (u), (-cB) ■ z < 0} where B,B f are as m 
(|) and (U) and r G A c . By Theorem p78| , both Q u and Q T U are polytopes. Notice that if 



n T (u) = tc t (u') for two distinct vectors u, u' G N n , then = Q T U ,. 

Lemma 3.4. (i) A lattice point u is in O c if and only if Q u D r L n ~ d = {0}. 

(ii) If u G O c , then the group relaxation G T (Au) solves the integer program IPa, c {Au) if and 

only ifQ f u C]Z n - d = {0}. 

Proof, (i) The lattice point u belongs to O c if and only if u is the optimal solution to 
IPa, c {Au) which is equivalent to G Z n ~ d being the optimal solution to the reformulation 
(|) of IPa, c (Au). Since c is generic, the last statement is equivalent to Q u PI Z n_d = {0}. 
The second statement follows from (i) and the fact that @ solves (f|) if and only if they 
have the same optimal solution. □ 

In order to state the coming results, it is convenient to assume that the vector u in (|j) 
and (Bp is the optimal solution to IPA, c {b). For an element u G O c and a face r of A c let 
S(u, t) be the affine semigroup w + N(ej : i G r) C N n where e» denotes the ith unit vector of 
W 1 . Note that S(u, r) is not a semigroup if u 7^ (since G" S^w, r)), but is a translation of 
the semigroup N(ej : i G r). We use the adjective affine here as in an affine subspace which 
is not a subspace but the translation of one. Note that if v G S(u, r), then ir T (v ) = 7t t (m). 

Lemma 3.5. For w G O c and a /ace r of A c , the affine semigroup S(u,t) is contained in 
O c if and only if G T (Au) solves IPa, c {Au). 



Proof. Suppose S(u, r) C O c . Then by Lemma [O] (i), for all v G S(u, r), 

Q v = {ze W l - d : B T z < 7T f (v), B f z < tt t («), (-c5) ■ 2 < 0} n Z n - d = {0}. 

Since TT f (v) can be any vector in Nl T l, Q f u n Z"- d = {0}. Hence, by Lemma U (ii), G T (Au) 
solves IPa, c {Au). 
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If v G S(u,t), then n T (u) = n T (v), and hence Q T U = Q T V . Therefore, if G T (Au) solves 
IP Ac (Au), then {0} = Q f u n Z n ~ d = Q f v n Z n - rf for all u G r). Since is a relaxation 
of Q„, n Z n ' d = {0} for all v G S(m, r) and hence by Lemma (i), r) C £> c . □ 

Lemma 3.6. Forw 6 C and a face r of A c , G T (Au) solves IPa, c {Au) if and only ifG T (Av) 
solves IPa, c (Av) for all v G S(u, r). 

Proof, live S(u, t) and G T (Au) solves IP A , c (Au), then as seen before, {0} = Ql n Z n_d = 
Qj;nZ n - d for allu G S(w,r). By Lemma (ii), G r (Aw) solves JP AjC (At;) for all u G S(u,t). 
The converse holds for the trivial reason that u G S(u, r). □ 

Corollary 3.7. For u G O c and a face t of A c , the affine semigroup S(u,r) is contained in 
O c if and only if G T (Av) solves IPa iC (Av) for all v G S(u, r). 

Since n T (u) determines the polytope Q T U = Ql for all v G S(u,r), we could have assumed 
that supp(u) C f in Lemmas |3.5| and |3.6| . 

Definition 3.8. For r G A c and u G O c , (u, r) is called an admissible pair of O c if 

(i) the support of u is contained in f , and 

(ii) S(u,t) C O c or equivalently, G T (Av) solves IPa, c {Av) for all v G S(u,t). 

An admissible pair (w, r) is a standard pair of O c if the affine semigroup S(u, r) is not prop- 
erly contained in S(v, r') where (v, r') is another admissible pair of O c . 



Example |3.2| continued. From Figure ^|, one can see that the standard pairs of O c 

((1,O,O),0) ((1,3,O),0) ((0,0,0), {3}) 

((2,O,O),0) ((2,3,O),0) ((0,1,0), {3}) 

((3,0,0), 0) ((1,4,0), 0) ((0,2,0), {3}) 

((1,1,O),0) ((1,5,O),0) ((0,3,0), {3}) 

((2,1,O),0) and ((0,4,0), {3}) 

((3,1,O),0) ((0,5,0), {3}) 

((1,2,O),0) ((0,6,0), {3}) 

((2,2,O),0) ((0,7,0), {3}) 



arc 



□ 



Definition 3.9. For a face r of A c and a lattice point u G N n , we say that the polytope Q T U 
is a standard polytope of IPa, c if Q T U H Z n ~ d = {0} and every relaxation of Q T U obtained by 
removing an inequality in B T z < tt t (u) contains a non-zero lattice point. 

Figure [| is a diagram of a standard polytope Q T U . The dashed line is the boundary of the 
half space (—cB) ■ z < while the other lines are the boundaries of the half spaces given by 
the inequalities in B T z < tc t (u). The origin is the only lattice point in the polytope and if 
any inequality in B T z < n T (u) is removed, a lattice point will enter the relaxation. 

We re-emphasize that if Q T U is a standard polytope, then Q T U , is the same standard polytope 
if 7iv (m) — 7r r (u'). Hence the same standard polytope can be indexed by infinitely many 
u G N n . We now state the main result of this section which characterizes associated sets in 
terms of standard pairs and standard polytopes. 
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Figure 5. A standard polytope. 



Theorem 3.10. The following statements are equivalent: 

(i) The admissible pair (u, r) is a standard pair of O c . 

(ii) The polytope Q T U is a standard polytope of IPa, c - 
(Hi) The face r of A c is associated to IPa,c- 

Proof, (i) (ii): The admissible pair (u,r) is standard if and only if for every i G r, 
there exists some positive integer TOj and a vector v G S(u, r) such that v + G Af c . (If 
this condition did not hold for some i G r, then (u 1 , r U {i}) would be an admissible pair 
of O c such that S(u', r U {i}) contains S(u, r) where u' is obtained from u by setting the 
ith component of u to zero. Conversely, if the condition holds for an admissible pair then 
the pair is standard.) Equivalently, for each i G r, there exists a positive integer rrii and a 
v G S(u, r) such that Q T v+m . e . = Q T u+miei contains at least two lattice points. In other words, 
the removal of the inequality indexed by % from the inequalities in B T z < ir T (u) will bring 
an extra lattice point into the corresponding relaxation of Q T U . This is equivalent to saying 
that Q T U is a standard polytope of IPa, c - 



(i) (Hi): Suppose (u,t) is a standard pair of O c . Then S(u,r) C O c and G T (Au) solves 
IPa, c (Au) by Lemma Suppose G T (Am) solves JPa,c(^4w) for some face r' G A c such 
that t C t' . Lemma |3]^ then implies that S(u, t') lies in O c . This contradicts the fact that 
(u, t) was a standard pair of O c since S(u, t) is properly contained in S(u, r') corresponding 
to the admissible pair (u, t') where u is obtained from u by setting Ui = for all i G r'\r. 

To prove the converse, suppose r is associated to IPa,c- Then there exists some b G PL4 
such that G T (b) solves IPA, c (b) but G T ' (b) does not for all faces r' of A c containing r. Let 
u be the unique optimal solution of IPA, c (b). By Lemma |3.5| , S(u,t) C (9 c . Let m G N" 
be obtained from u by setting Mj = for all i G r. Then (^(Afi) solves JP J 4,c(A^) since 
= Q T u- Hence S(u,r) C (9 C and (m, r) is an admissible pair of O c . Suppose there exists 
another admissible pair (w,a) such that S(u,t) C ^(w,^). Then r C a. If r = <r then 
S'('U, r) and 5 , (w, a) are both orthogonal translates of N(ej : i G r) and hence S(u, r) 
cannot be properly contained in S(w, a). Therefore, r is a proper subset of a which implies 
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that S(u,a) C O c . Then, by Lemma |3~5|, G a (Au) solves IPa >c (Au) which contradicts that r 
was an associated set of I Pa c - D 



Example |3.2| continued. In Example |3.2| we can choose B to be the 3x2 matrix 

-1 



B 



4 




1 



The standard polytope defined by the standard pair ((1, 0, 0), 0) is hence 



{(zi,z 2 ) e 



-zi + Az 2 < 1, 2zi < 0, -z x - z 2 < 0, 9801^1 - 40001z 2 < 0} 



while the standard polytope defined by the standard pair ((0, 2, 0), {3}) is 

{{z x , z 2 ) G M 2 : -zi + 4^ < 0, 2z x < 2, 9801*1 - 40001z 2 < 0}. 

The associated sets of IPa.c in this example are and {3}. There are twelve quadrangular 
and eight triangular standard polytopes for this family of knapsack problems. □ 

Standard polytopes were introduced in |22| , and the equivalence of parts (i) and (ii) of 
Theorem |3.10| was proved in |22|, Theorem 2.5]. Under the linear map 4>a '■ N n —>■ NA where 
u i — > Au, the affine semigroup S(u, r) where (u, r) is a standard pair of O c maps to the 
affine semigroup Au + NA T in NA. Since every integer program in IPa, c is solved by one of 
its group relaxations, O c is covered by the affine semigroups corresponding to its standard 
pairs. We call this cover and its image in NA under (pA the standard pair decompositions 
of O c and NA, respectively. Since standard pairs of O c are determined by the standard 
polytopes of IPa,ci the standard pair decomposition of O c is unique. The terminology used 
above has its origins in |58J which introduced the standard pair decomposition of a monomial 
ideal. The specialization to integer programming appear in [^3J and §12. D]. The 
following theorem shows how the standard pair decomposition of O c dictates which group 
relaxations solve which programs in IPa,c- 

Theorem 3.11. Let v be the optimal solution of the integer program JP^ C (6). Then the 
group relaxation G T (Av) solves IPa, c {Av) if and only if there is some standard pair (u,t') 
of O c with t C t' such that v belongs to the affine semigroup S(u, r'). 

Proof. Suppose v lies in S(u, r') corresponding to the standard pair (u, r') of O c . Then 
S{v,t') C O c which implies that G T (Av) solves IPa, c {Av) by Lemma p.5| . Hence G T (Av) 
also solves IPa, c {Av) for all r C r'. 

To prove the converse, suppose r' is a maximal element in the subcomplex of all faces r 
of A c such that G T (Av) solves IPa, c (Av). Then r' is an associated set of IPa,c- In the proof 
i) in Theorem p.!0| , we showed that (-0, r') is a standard pair of O c where v is 
for all % G r' . Then v G S(v, r'). □ 



of (Hi) - 
obtained from v by setting v. 



Example |3.2| continued. The eight standard pairs of O c of the form (*, {3}), map to the 
eight affine semigroups: 

N[8], (5 + N[8]), (10 + N[8]), (15+N[8]), (20 + N[8]), (25 + N[8]), (30 + N[8]) and (35 + N[8]) 
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contained in NA = N [2, 5, 8] C N. For all right hand side vectors b in the union of these 
sets, the integer program IP^ tC (b) can be solved by the group relaxation G^(b). The twelve 
standard pairs of the from (*, 0) map to the remaining finitely many points 

2, 4, 6, 7, 9, 11, 12, 14, 17, 19, 22 and 27 

of N [2, 5, 8]. If b is one of these points, then JP J 4 jC (6) can only be solved as the full integer 
program. In this example, the regular triangulation A c = {{3}}. Hence G^(b) is a Gomory 
relaxation of IP^ c (b). □ 



For most b G NA, the program IP^ c (b) is solved by one of its Gomory relaxations, or 



equivalently, by Theorem 3.11 , the optimal solution v of IPa, c {P) lies in S(*, a) for some 
standard pair (*, a) where a is a maximal face of A c . For mathematical versions of this in- 
formal statement, see 0, Proposition 12.16] and |13|, Theorems 1 and 2]. Roughly speaking, 
these right hand sides are away from the boundary of cone(A). (This was seen in Exam- 
ple |3.2j above, where for all but twelve right hand sides, IPA, c {b) was solvable by the Gomory 
relaxation G^(b). Further, these right hand sides were toward the boundary of cone(A), 
the origin in this one-dimensional case.) For the remaining right hand sides, IPA, c {b) can 
only be solved by G T (b) where r is a lower dimensional face of A c - possibly even the empty 
face. An important contribution of the algebraic approach here is the identification of the 
minimal set of group relaxations needed to solve all programs in the family IPa, c and of the 
particular relaxations necessary to solve any given program in the family. 



4. Arithmetic Degree 

For an associated set r of IPa, c there are only finitely many standard pairs of O c indexed 
by r since there are only finitely many standard polytopes of the form Q T U . Borrowing 
terminology from |38|| , we call the number of standard pairs of the form (•, r) the multiplicity 
of r in O c (abbreviated as mult{r)). The total number of standard pairs of O c is called 
the arithmetic degree of O c . Our main goal in this section is to provide bounds for these 
invariants of the the family IPa.c and discuss their relevance. We will need the following 
interpretation from Section 3. 

Corollary 4.1. The multiplicity of the face r of A c in O c is the number of distinct standard 
polytopes of IPa, c indexed by f, and the arithmetic degree of O c is the total number of 
standard polytopes of IPa, c - 

Proof. This result follows from Theorem [3.10| . □ 

Example |3.2| continued. The multiplicity of the associated set {3} is eight while the 
empty set has multiplicity twelve. The arithmetic degree of O c is hence twenty. □ 



If the standard pair decomposition of O c is known, then we can solve all programs in IPa, c 
by solving (arithmetic degree)-many linear systems as follows. For a given b G NA and a 
standard pair (u,t), consider the linear system 



(6) 



Af7r T (u) + A T x = b, or equivalently, A T x = b — A f ir T (u). 



AN ALGEBRAIC PERSPECTIVE OF GROUP RELAXATIONS 



17 



As r is a face of A c , this linear system can be solved uniquely for x. Since the optimal solution 
of IPa,c(1>) hes in S(w, a) for some standard pair (w, a) of O c , at least one non-negative and 
integral solution for x will be found as we solve the linear systems (^j) obtained by varying 
(u, t) over all the standard pairs of O c . If the standard pair (u, r) yields such a solution 
v, then (tt t (u),v) is the optimal solution of IPA <c {b). This pre-processing of IPa, c has the 
same flavor as ||27f . The main result in []27| is that given a coefficient matrix A G M mxn and 
cost vector c, there exists floor functions fi, ■ ■ ■ , fk '■ K m — » ^ n such that for a right hand 
side vector b, the optimal solution of the corresponding integer program is the one among 
/i (&),... , fkip) that is feasible and attains the best objective function value. The crucial 
point is that this algorithm runs in time bounded above by a polynomial in the length of 
the data for fixed n and j, where j is the affine dimension of the space of right hand sides. 
Given this result, it is interesting to bound arithmetic degree. 

The second equation in @ suggests that one could think of the first arguments u in 
the standard pairs (u, r) of O c as "correction vectors" that need to be applied to find the 
optimal solutions of programs in IPa, c - Thus the arithmetic degree of O c is the total number 
of correction vectors that are needed to solve all programs in IPa, c - The multiplicities of 
associated sets give a finer count of these correction vectors, organized by faces of A c . If the 
optimal solution of IPA, c {b) lies in the affine semigroup S(w, a) given by the standard pair 
(w, a) of O c , then w is a correction vector for this b as well as all other 6's in (Aw + N^). 
One obtains all correction vectors for IPa, c by solving the (arithmetic degree)-many integer 
programs with right hand sides Au for all standard pairs (u, r) of O c . See [[PJ for a similar 
result from the classical theory of group relaxations. 

In Example |3]^, A c = {{3}} and both its faces {3} and are associated to IPa, c - In 
general, not all faces of A c need be associated sets of IPa, c and the poset of associated sets 
can be quite complicated. (We will study this poset in Section 5.) Hence, for r G A c , 
mult(r) = unless r is an associated set of IPa,c- We will now prove that all maximal faces 
of A c are associated sets of IPa,c- Further, if a is a maximal face of A c then mult(o~) is 
the absolute value of det(A a ) divided by the g.c.d. of the maximal minors of A. This g.c.d 
is non-zero since A has full row rank. If the columns of A span an affine hyperplane, then 
the absolute value of det(A a ) divided by the g.c.d. of the maximal minors of A is called the 
normalized volume of the face a in A c . We first give a non-trivial example. 

Example 4.2. Consider the rank three matrix 



and the generic cost vector c = (21, 6, 1, 0, 0, 0). The first three columns of A generate 
cone(A) which is simplicial. The regular triangulation 



is shown in Figure |^ as a triangulation of conv(A). The six columns of A have been labeled 
by their column indices. The arithmetic degree of O c in this example is 70. The following 
table shows all the standard pairs organized by associated sets and the multiplicity of each 
associated set. Note that all maximal faces of A c are associated to IPa, c - The g.c.d. of the 



.4 



5 
5 
5 



2 
1 

2 



1 
4 2 
3 



A c = {{1,3,4}, {1,4,5}, {2,5,6}, {3,4,6}, {4,5,6}} 
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FIGURE 6. The regular triangulation A c for Example [4.2| . 

maximal minors of A is five. Check that mult(a) is the normalized volume of a whenever a 
is a maximal face of A c . 



T 


standard pairs (-,r) 


mult(r) 


11 3 4-1 






11 4 51 


(3e 2 , •), (2e 2 + e 3 , •) 


8 


{2,5,6} 


(0,-), (e 3 ,-), (2e 3 ,-) 


3 


{3,4,6} 


(0,-), (e B ,-), (2e 5 ,-), (3e 5 ,-) 


4 


{4,5,6} 


(0,-), (e 3 ,-), (2e 3 ,-), (3e 3 ,-), (4e 3 , •) 


5 


{1,4} 


(e 3 + 2e 5 + e 6 , •), (2e 3 + 2e 5 + e 6 , •), (2e 3 + 2e 5 , •), 
(2e 3 + 3e 5 ,-), (2e 3 + 4e 5 , •) 


5 


{1,5} 


(e 2 + e 6 , •), (2e 2 + e 6 , ■), (3e 2 + e 6 , •) 


3 


{2,5} 


(e 3 + e 4 , ■), (e 4 , ■), (2e 4 , •) 


3 


{3,4} 


(e 2 , •), (ei + e 2 , ■), (ei + 2e 5 , •), (ei + 2e 5 + e 6 , ■), (e 2 + e 6 , •), 


5 


{3,6} 


(e 2 , ■), (e2 + e 5 , •) 


2 


{4,5} 


(e 2 + 2e 3 , •), (e 2 + 3e 3 , ■), (2e 2 + 2e 3 , ■), (3e 2 + e 3 , ■), (4e 2 , •) 


5 


{5,6} 


(e 2 + 3e 3 , •) 


1 


{1} 


(e 2 + e 3 + e 6 , •), (e 2 + e 3 + e 5 + e 6 , •), (e 2 + 2e 6 , •), 
(e 2 + e 3 + 2e 6 , •), (2e 2 + 2e 6 , •), (e 2 + e 3 + 2e 5 + e 6 , •) 


6 


{3} 


(ei + e 2 + e 6 , •), (ei + e 2 + 2e 6 , •) 


2 


{4} 


(ei + e 2 + 2e 3 + e 5 , •), (ei + e 2 + 2e 3 + 2e 5 , •), 
(ei + e 2 + 2e 3 + 3e 5 , •), (ei + e 2 + 2e 3 + 4e 5 , •), 
(ei + 3e 3 + 3e 5 , •), (d + 3e 3 + 4e 5 , •) 


6 


{0} 


(ei + e 2 + 2e 3 + e 5 + e 6 , •), (ei + e 2 + 2e 3 + 2e 5 + e 6 , •), 
(ei + 2e 2 + e 3 + e 6 , •), (e x + 2e 2 + e 3 + e 5 + e 6 , •), 
(ei + 2e 2 + e 3 + 2e 5 + e 6 , •), (e a + 2e 2 + e 3 + 2e 6 , •), 
(ei + 3e 2 + 2e 6 ,-) 


7 




Arithmetic Degree 
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Observe that the integer program IP^ c (b) where b = A(ei + e2 + e^) is solved by G T (b) with 
r = {1, 4, 5}. By Proposition p73| , Gomory's relaxation of IPa,c(P) is indexed by o = {4, 5, 6} 
since b lies in the interior of the face cone{A a ) of A c . However, neither this relaxation nor 
any nontrivial extended relaxation solves IPA, c {b) since the optimal solution e\ + C2 + e% is 
not covered by any standard pair (*, r) where r is a non-empty subset of {4, 5, 6}. □ 



Theorem 4.3. For a set a C {1, . . . , n}, (0, er) is a standard pair of O c if and only if a is 
a maximal face of A c . 

Proof. If a is a maximal face of A c , then by Definition |2.1| , there exists y G M. d such that 
yA a = c a and yA„ < c 5 . Thenp = c 5 — yA 5 > a,iadpB a = (c^—yA d .)B a = c 5 B a ' —yA 5 B a = 
c^B a + yA a B a = c^B a + c a B a = cB. Hence there is a positive dependence relation among 
(—cB) and the rows of B° '. Since a is a maximal face of A c , |G?et( J 4 cr ) | ^ 0. However, 
\det(B a ) | = |c?et(y4 cr ) | which implies that \det(B cr )\ ^ 0. Therefore, {—cB) and the rows of 
B° span W 1 - 6 - positively This implies that Q% = {z G R n ~ d : B 5 "z < 0, (-cB) ■ z < 0} is 
a polytope consisting of just the origin. If any inequality defining this simplex is dropped, 
the resulting relaxation is unbounded as only n — d inequalities would remain. Hence Qq is 
a standard polytope of IPa, c and by Theorem |3.10| , (0, a) is a standard pair of O c 



Conversely, if (0, a) is a standard pair of O c then Q% is a standard polytope of IPa,c- 
Since every inequality in the definition of Qq gives a halfspace containing the origin and Qq 
is a polytope, Qq = {0}. Hence there is a positive linear dependence relation among (—cB) 
and the rows of B° ' . If |cr| > n — d, then Qq would coincide with the relaxation obtained 
by dropping some inequality from those in B° ' z < 0. This would contradict that Qq was a 
standard polytope and hence |er| = d and a is a maximal face of A c . □ 

Corollary 4.4. Every maximal face of A c is an associated set of IPa, c - 



For Theorem |4.5| and Corollary [05] below we assume that the g.c.d. of the maximal minors 
of A is one which implies that *LA = Z, d . 

Theorem 4.5. If a is a maximal face of A c then the multiplicity of a in O c is \det(A a )\. 

Proof. Consider the full dimensional lattice C a = 7r CT (£) = {B° ~z : z G Z n ~ d } in Z n ~ d . 
Since the g.c.d. of the maximal minors of A is assumed to be one, the lattice L a has index 
\det(B a ) \ = |c?et(A CT )| in Z n ~ d . Since L a is full dimensional, it has a strictly positive element 
which guarantees that each equivalence class of Z n ~ d modulo C a has a non-negative member. 
This implies that there are |<iet(A (T )| distinct equivalence classes of W"~ d modulo C a . Recall 
that if u is a feasible solution to JP J 4 jC (6) then 

G a {b) = minimize {c& ■ x„ : x„ = u„ {mod C a ), x„ G N n ~ d }. 

Since there are |c?et(y4 cr ) | equivalence classes of M n ~ d modulo C a , there are ((ie^Ao-)! distinct 
group relaxations indexed by a. The optimal solution of each program becomes the right 
hand side vector of a standard polytope (simplex) of IPa, c indexed by o. Since no two 
optimal solutions are the same (as they come from different equivalence classes of W L ~ d 
modulo C a ), there are precisely |<iet (y4 CT ) | standard polytopes of IPa, c indexed by a. □ 
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Corollary 4.6. The arithmetic degree of O c is bounded below by the sum of the absolute 
values of det(A a ) as a varies among the maximal faces of A c . 

A primary ideal J in k[xi, . . . , x n ] is a proper ideal such that fg G J implies either / G J 
or g l G J for some positive integer t. A prime ideal J of k[xi, . . . ,x n ] is a proper ideal 
such that fgEJ implies that either / G J or g G J. A primary decomposition of an 
ideal J in . . . ,x n ] is an expression of / as a finite intersection of primary ideals in 



k[xi, . . . , x n ]. Lemma 3.3 in [38| shows that every monomial ideal M in k[x\, . . . , x n ) admits 
a primary decomposition into irreducible primary ideals that are indexed by the standard 
pairs of M. The radical of an ideal I C k[xi, ... , x n ] is the ideal a/7 := {/ 6 S : /* G 
J, for some positive integer t}. Radicals of primary ideals are prime. The radicals of the 
primary ideals in a minimal primary decomposition of an ideal I are called the associated 
primes of /. This list of prime ideals is independent of the primary decomposition of the 
ideal. The minimal elements among the associated primes of / are called the minimal 
primes of / while the others are called the embedded primes of /. The minimal primes of 
/ are precisely the defining ideals of the isolated components of the zero-set or variety of / 
while the embedded primes cut out embedded subvarieties in the isolated components. See 



a textbook in commutative algebra like |H| for more details. 

A face r of A c is an associated set of IPa, c if and only if the monomial prime ideal 
p T := (xj : j ^ t) is an associated prime of the ideal in c (lA). Further, p a is a minimal prime 
of in c {J.A) if and only if a is a maximal face of A c . Hence the lower dimensional associated 
sets of IPa, c index the embedded primes of in c {lA)- The standard pair decomposition of a 
monomial ideal was introduced in JJS] to study its associated primes. The multiplicity of 
an associated prime p T of in c (lA) is an algebraic invariant of in c (lA), and shows that 
this is exactly the number of standard pairs indexed by r. Similarly, the arithmetic degree 
of in c (lA) is a refinement of the geometric notion of degree and shows that this number 
is the total number of standard pairs of in c (lA)- These connections explain our choice of 
terminology. Theorem (Ol is a translation of the specialization of Lemma 3.5 in PSJ to 



toric initial ideals. We refer the interested reader to BO, §8 and §12. D] and PHI §3] for the 



algebraic connections. Theorem [15] is a staple result of toric geometry and also follows from 



fl3| , Theorem 1]. It is proved via the algebraic technique of localization in [3£| Theorem 8.8]. 

Theorem [4.5| gives a precise bound on the multiplicity of a maximal associated set of 
IPa,c, which in turn provides a lower bound for the arithmetic degree of O c in Corollary [4.6| . 
No exact result like Theorem [L5] is known when r is a lower dimensional associated set of 
IPa,c- Such bounds would provide a bound for the arithmetic degree of O c . Bounds on 
the arithmetic degree of a general monomial ideal in terms of its dimension and minimal 



generators can be found in |}8, Theorem 3.1]. One hopes that stronger bounds are possible 
for toric initial ideals. We close with a first attempt at bounding the arithmetic degree of 
O c (under certain non-degeneracy assumptions). This result is due to Ravi Kannan, and its 



simple arguments are along the lines of proofs in [26[] and |28j . 

Suppose S G Z mxn and u G N m are fixed and K u := {x G IR n : Sx < u] is such that 
K u fl Z™ = {0} and the removal of any inequality defining K u will bring in a non-zero lattice 
point into the relaxation. Let denote the ith row of S, M := raax||s^||i and Afc(S') and 
5k{S) be the maximum and minimum absolute values of the k x k sub determinants of 5*. We 
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will assume that S n (S) 7^ which is a non-degeneracy condition on the data. We assume 



this set up in Theorem [4.8| and Lemmas 4.9 and 4.10 



Definition 4.7. If K is a convex set and v a non-zero vector in R n , the width of K along 

v, denoted as width v (K) is max {v ■ x : 16 K} — min {v ■ x : 16 K}. 



Note that width v (K) is invariant under translations of K. 
Theorem 4.8. If K u is as above then < m < 2M(n + 2)^^. 

Lemma 4.9. If K u is as above then for some t, 1 < t < m, width s ( t )(K u ) < Min + 2). 



Proof. Clearly, K u is bounded since otherwise there would be a non-zero lattice point on an 
unbounded edge of K u due to the integrality of all data. Suppose width 8 (t)(K u ) > M{n + 2) 
for all rows s^> of S. Let p be the center of gravity of K u . Then by a property of the center 
of gravity, for any x G K u , (l/(n + l))th of the vector from p to the reflection of x about p 
is also in K u , i.e., (1 + -^i)p — ^p[ x K u . Fix i, 1 < i < m and let xq minimize ■ x over 
K u . By the definition of width, we then have — ■ Xo > M(n + 2) which implies that 

(7) s {i) ■ x <Ui- M(n + 2). 

Now s w ■ ((1 + ^i)p - ^1^0) < ^ implies that 

/^v m ,n + l s sW ■ xr, 

8 s W -p<^ (— - +— 

n + 2 n + 2 

Combining (0) and (Eft we get 

(9) ■ p < Ui — M. 

Let g = [pj be the vector obtained by rounding down all components of p. Then p = q + r 
where < < 1 for all j = 1, . . . , n, and by (Q), • (g + r) < u\ — M which leads to 
S W . g _|_ ( s ( l ) . r 4- ]Vf) < Mi . Since M = max\\s^ | |i, 

(10) — M < s {i) -r<M. 

and hence, • g < Repeating this argument for all rows of S, we get that q G -ft'u. 
Similarly, if q' = \p] is the vector obtained by rounding up all components of p, then p = q'—r 
where < r j < 1 for all j — 1, . . . , n. Then (^j) implies that s 1 -^ • (q' — r) < Ui — M which 
leads to s^' ■ q' + (M — s^' ■ r) < u^. Again by (PT|), ■ q' < and hence q' G K u . Since 
q 7^ q', at least one of them is non-zero which contradicts that K u fl Z™ = {0}. □ 

Lemma 4.10. For any two rows s^\s^ of S, width s (€,{K u ) < 2 width s u) (K u ) . 

Proof. Without loss of generality we may assume that j = n + 1. Since K u is bounded, 
width s ( 3 )(K u ) is finite. Suppose the minimum of s^> ■ x over K u is attained at v. Since 
translations leave the quantities in the lemma invariant, we may prove the lemma for the 
body K u i obtained by translating K u by —v. Now s^> -x is minimized over K u i at the origin. 
By LP duality, there are n linearly independent constraints among the m defining K u i such 
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that the minimum of s^ n+1 ^ ■ x subject to just these n constraints is attained at 0. After 
renumbering the inequalities if necessary, assume these n constraints are the first n. Let 

D = {x : s {l) ■ x < u[, I = 1, 2, . . . , n + 1} 

where of course u[ = u' 2 = ■ ■ ■ = u' n = 0. Then by the above, D is a bounded simplex. 
Since D contains K u i, it suffices to show that for each i, 

(11) widthMD) < 2(^-)width s(n+1) (K u ,) = 2(^§K +1 . 

We show that for each vertex q of D, |sW . g| < ) u ' n +i which will prove flTT|). This is 

clearly true for q = 0. Without loss of generality assume that vertex q satisfies s"'-q = u\ 
for I — 2, 3, . . . , n + 1. Since the determinant of the submatrix of 5* consisting of the rows 
s^ 2 -*, . . . , s( n+1 ) is not zero, for any i there exists rationals A; such that = Ya=2 ^i s ^ ■ By 
Cramer's rule, |A;| < ( ) • Therefore, sW • q = J2™=2 ^i- 3 ^ ' Q = ^ lU 'i = ^n+iu' n+1 

since u[ = for I = 2, . . . ,n. This proves that 

□ 



k w • g| = = |A n+ i|< +1 < ( \y o ( )u' n+v 



Proof of Theorem \^.8j . From Lemmas and [4.10| it follows that for any i, 1 < i < m, 



width s( i)(K u ) < 2(j^§j)M(n + 2) = 2M(n + 2)(^^). Since E K u , min{s® ■ x : x E 
Ku] < while max{s^ ■ x : x E K u } = U{. Therefore, Ui = Ui — < width s (i)(K u ) and 
hence, < u { < 2M(n + 2)(^^) for all 1 < i < m. □ 



Reverting back to our set up, let B 



B 

-cB 



Suppose K u is the standard polytope 



Q f u . By Theorem p|, < u { < 2M{n -d+ 2)(j^§± 



Corollary 4.11. If no maximal minor of B is zero, then the arithmetic degree of O c is at 
most (2M(n -d + 2)(^=^)V. 

The above arguments do not use the condition that the removal of an inequality from K u 
will bring in a lattice point into the relaxation. Further, the bound is independent of the 
number of facets of K u , and Corollary |4.11| is straightforward. Thus, further improvements 



may be possible with more effort. However, apart from providing a bound for arithmetic 
degree, these proofs have the nice feature that they build a bridge to techniques from the 
geometry of numbers that have played a central role in theoretical integer programming as 
seen in the work of Kannan, Lenstra, Lovasz, Scarf and others. See [29|] for a survey. 



5. The Chain Theorem 

We now examine the structure of the poset of associated sets of IPa, c which we denote as 
Assets(IPA, c )- All elements of Assets (IPa c ) are faces of the regular triangulation A c and 
the partial order is set inclusion. Theorem PO provides a first result. 
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Figure 7. Lower dimensional associated sets of Example |4.2| except the 
empty set. 



Corollary 5.1. The maximal elements of Assets(IPA tC ) are the maximal faces of A c 



Example |4.2| continued. The lower dimensional associated sets of this example (except 



the empty set) are the thick faces of A c shown in Figure [7]. □ 

Despite the seemingly chaotic structure of Assets(IPA, c ) beyond its maximal elements, it 
has an important structural property that we now explain. 

Theorem 5.2. [The Chain Theorem] If t G A c is an associated set of IPa,c and \t\ < d 
then there exists a face r' G A c that is also an associated set of IPa, c with the property that 
t C t' and \t'\t\ = 1. 

Proof. Since r is an associated set of IPa,c, by Theorem |3.10| , O c has a standard pair of the 
form (v,t) and Q T V = {z G W l ^ d : B T z < tt t (v), {—cB) ■ z < 0} is a standard polytope of 
IPa, c - Since |r| < d, r is not a maximal face of A c and hence by Theorem |4.3| , v ^ 0. For 
each i G r, let R l be the relaxation of Q T V obtained by removing the ith inequality bi ■ z < v { 
from B T z < tc t (v), i.e., 

R} := {z G W- d : B f \ {i} z < TT TU{i} (v), (-cB) -z<0}. 

Let E % := R l \Ql. Clearly, E l f~)Ql = 0, and, since the removal of b^-z < Vi introduces at least 
one lattice point into R x , E l r\Z n ~ d ^ 0. Let z* be the optimal solution to minimize {(—cB) ■ 
z : z G E 1 H Z n ~ d } if the program is bounded. This integer program is always feasible since 
E % fl 7L n ~ d 7^ 0, but it may not have a finite optimal value. However, there exists at least 
one i G r for which the above integer program is bounded. To see this, pick a maximal 
simplex a G A c such that r C a. The polytope {z G R n ~ d : B°z < 7T a (v), (-cB) ■ z < 0} 
is a simplex and hence bounded. This polytope contains all E l for i G c\t, and hence all 
these E l are bounded and have finite optima with respect to (— cB) ■ z. We may assume that 
the inequalities in B T z < 7i T (v) are labeled so that the finite optimal values are ordered as 
(-cB) ■ z\ > (-cB) ■ z* 2 > ■ ■ ■ > (-cB) ■ z* where {1, 2, . . . ,p} C f . 
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Claim: Let N 1 := {z G R n ~ d : B f \^z < n TU{1} (v), (-cB) ■ z < (-cB) ■ z\}. Then z\ is 
the unique lattice 'point in N 1 and the removal of any inequality from B T ^^z < 7r rU {i}(f) 
will bring in a new lattice point into the relaxation. 

Proof. Since z\ lies in R 1 , = (—cB) ■ > (—cB) ■ z\. However, > {—cB) ■ z\ since 
otherwise, both z\ and would be optimal solutions to minimize{(—cB) ■ z : z G -R 1 } con- 
tradicting that c is generic. Therefore, 
N 1 = R 1 n {z G R n ~ d : (-cB) ■ z < (-cB) ■ z{} 

= (E 1 U Ql) C]{zE R n ~ d : (-cB) ■ z < (-cB) ■ z{} 

= (E 1 n {z G R n ~ d : (-cB) ■ z < {-cB) ■ z{}) 
v U (Ql n{ze R n ' d : {-cB) ■ z < (-cB) ■ z{}). 
Since c is generic, z\ is the unique lattice point in the first polytope and the second polytope 
is free of lattice points. Hence z\ is the unique lattice point in iV 1 . The relaxation of N 1 got 
by removing bj ■ z < Vj is the polyhedron iV 1 U (E j fl {z G M n ~ d : (—cB) ■ z < (—cB) ■ z{}) 
for j G f and j ^ 1. Either this is unbounded, in which case there is a lattice point z in this 
relaxation such that (—cB) ■ z\ > (—cB) ■ z, or (if j < p) we have (—cB) ■ z\ > (—cB) ■ z* 
and z* lies in this relaxation. <) 



Translating A^ 1 by -z\ we get Q T V ) { } := {z G R n - d : (-cB) ■ z < 0, B^h < v'} where 
t>' = tt t ij^(v) — B T \^zl > since z\ is feasible for all inequalities except the first one. Now 
Q f v } {1} n % n ~ d = {0}, and hence (V, r U {!}) is a standard pair of O c . □ 



Example |4.2| continued. The empty set is associated to IPa, c an d C {1} C {1,4} C 



{1,4,5} is a saturated chain in Assets (IPa, c ) that starts at the empty set. □ 
In algebraic language, the chain theorem says that the associated primes of in c (Ij{) occur 



in saturated chains. This was proved in p2| , Theorem 3.1]. When the cost vector c is not 
generic, in c (lA) is no longer a monomial ideal, and its associated primes need not come in 
saturated chains. See |22], Remark 3.3] for such an example. An important open question 



in the algebraic study of integer programming is to characterize all monomial ideals that 
can appear as the initial ideal (with respect to some generic cost vector) of a toric ideal. 
In our set up this amounts to characterizing all down sets in N n that can appear as the 
set of optimal solutions to a family IPa.c, where A and c satisfy the assumptions from 
Section 2. Theorem |5.2| imposes the necessary condition that the poset of sets indexing the 
standard pairs of the down set have the chain property. Unfortunately, this is not sufficient 
to characterize down sets of the form O c . See |$0| for another class of monomial ideals that 



also have the chain property. 

Since the elements of Assets(IPA, c ) are faces of A c , a maximal face of which is a <i-element 
set, the length of a maximal chain in Assets(IPA, c ) is at most d. We denote the maximal 
length of a chain in Assets(IPA, c ) by length(Assets(I Pa,c)) ■ When n — d (the corank of A) 
is small compared to d, length(Assets(I Pa,c)) has a stronger upper bound than d. We use 
the following result of Bell and Scarf to prove the bound. 



Theorem 5.3. |33], Corollary 16.5 a] Let Ax < b be a system of linear inequalities in n 
variables, and let c G R n . If max {c ■ x : Ax < b, x G Z n } is a finite number, then max 
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{c ■ x : Ax < b, x G Z n } = max {c ■ x : A'a; < b',x G Z™} /or some subsystem A'x < b' of 
Ax < b with at most 2 n — 1 inequalities. 

Theorem 5.4. T/je length of a maximal chain in the poset of associated sets of IPa,c is at 
most min(d, 2 n ~ d — (n — d + 1)) . 

Proof. As seen earlier, length(Assets(I Pa,c)) < d. If v lies in O c , then the origin is the 
optimal solution to the integer program minimize {(—cB) ■ z : Bz < v, z G Z n ~ d }. By 
Theorem |5.3|, we need at most 2 n ~ d — 1 inequalities to describe the same integer program 



which means that we can remove at least n — (2 n — 1) inequalities from Bz < v without 
changing the optimum. Assuming that the inequalities removed are indexed by r, Q T V will be 
a standard polytope of IPa, c - Therefore, \t\ > n — (2 n ~ d — 1). This implies that the maximal 
length of a chain in Assets[lP A , c ) is at most d - (n - {2 n ' d - 1)) = 2 n ~ d -(n-d +1). □ 

Corollary 5.5. The cardinality of an associated set of IPa.c is at least max(0,n— (2 n ~ d — 1)). 

Corollary 5.6. If n — d = 2, then length(Assets(I Pa,c)) < 1- 

Proof. In this situation, 2 n ~ d - + = 4- (4 - 2 + 1) = 4- 3 = 1. □ 

We conclude this section with a family of examples for which length(Assets(IPA, c )) = 
2 n ~ d — (n — d+ 1). This is adapted from [p2|, Proposition 3.9] which was modeled on a family 
of examples from [|32]j . 

Proposition 5.7. For each m > 1, there is an integer matrix A of corank m and a cost 
vector cGZ" where n = 2 m — 1 such that length(Assets(I P A ,c)) = 2 m — (m + 1). 

Proof. Given m > 1, let B' = (&y) G ^( 2m - 1 ) xm be the matrix whose rows are all the 
{1, -l}-vectors in R m except v = (-1, -1, . . . , -1). Let B G Z ( - 2m+m -^ xm be obtained by 
stacking B' on top of — I m where I m is the m x m identity matrix. Set n = 2 m + m — 1, 
d = 2 m — 1 and A' = [I^B'] G Z d><n . By construction, the columns of B span the lattice 
{u G Z n : A'm = 0}. We may assume that the first row of B' is (1, 1, . . . , 1) G M. m . Adding 
this row to all other rows of A' we get A G N dxn with the same row space as A'. Hence the 
columns of B are also a basis for the lattice {u G Z™ : Au = 0}. Since the rows of B span 
Z m as a lattice, we can find a cost vector c G Z n such that (— cB) = v. 

For each row b{ of 5' set r« := 1(6^ : 6^ = 1}|, and let r be the vector of all r^s. By 
construction, the polytope Q := {z G IR m : B'z < r, —(cB) ■ z < 0} has no lattice points in 
its interior, and each of its 2 m facets has exactly one vertex of the unit cube in IR m in its 
relative interior. If we let Wi = rj — 1, then the polytope {z G M m : B'z < w, —(cB) ■ z < 0} 
is a standard polytope Q T U of IPa, c where r = {d+ 1, d + 2, . . . ,d + m = n} and w = n T (u). 
Since a maximal face of A c is a = (2 m — l)-element set and |r| = m, Theorem |57| implies 
that length(Assets(IPA. c )) > 2 m — 1 — m = 2 m — (m + 1). However, by Theorem [5.4| , 
length(Assets(IP A ,c)) = min(2 m - l,2 m - (m + 1)) = 2 m - (m + 1) since m > 1 by 
assumption. □ 
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Example 5.8. If we choose m = 3 then n = 2 m + m — 1 = 10 and d 



Constructing B' and A as in Proposition |5.7|, we get 



B' 
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1 ■ 




- 1 




















1 


1 


1 


-1 


1 


1 




1 


1 




















2 


2 


1 


-1 


1 




1 





1 














2 





2 


1 


1 


-1 
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The vector c = (11, 0, 0, 0, 0, 0, 0, 10, 10, 10) satisfies (-cB) = 
sets of IPa,c along with their multiplicities are given below. 



The associated 



T 


Multiplicity 


T 


Multiplicity 


{4,5,6, 7,8,9,10}* 


4 


{2,3,7,8,9,10} 


2 


{1,5,6,7,8,9,10} 


4 


{5,6,7,8,9,10}* 


1 


{3,4,6,7,8,9,10} 


4 


{4,5,6,7,8,9} 


1 


{2,3,4,6,7,9,10} 


2 


{2,4,7,8,9,10} 


2 


{2,3,4, 7,8,9,10} 


4 


{1,5,7,8,9,10} 


1 


{3,4,5,6,7,8,10} 


2 


{2,3,4,8,9,10} 


1 


{2,3,4,5,6,7,10} 


1 


{4,5,7,8,9,10} 


2 


{2,4,5,6,7,9,10} 


2 


{2,5,6,7,9,10} 


1 


{2,3,6,7,9,10} 


1 


{4,5,6,8,9,10} 


2 


{3,4,5,6,8,10} 


1 


{1,5,6,8,9,10} 


1 


{2,4,5,7,9,10} 


1 


{3,4,6,8,9,10} 


2 


{1,6,7,8,9,10} 


1 


{6,7,8,9,10}* 


1 


{3,5,6,7,8,10} 


1 


{7,8,9,10}* 


1 


{3,6,7,8,9,10} 


2 


{8,9,10}* 


1 



The elements in the unique maximal chain in Assets (I Pa 
length(Assets(I Pa, c )) = 2 3 — (3 + 1) = 4 as predicted by Proposition [T7 



are marked with a * and 

□ 



6. Gomory Integer Programs 



Recall from Definition 2~77 that a group relaxation G a (b) of IPA, c {b) is called a Gomory 
relaxation if o is a maximal face of A c . As discussed in Section 2, these relaxations are the 
easiest to solve among all relaxations of /Pa c (6). Hence it is natural to ask under what 
conditions on A and c would all programs in IP a c be solvable by Gomory relaxations. We 



study this question in this section. The majority of the results here are taken from [21 



Definition 6.1. The family of integer programs IPa, c is a Gomory family if, for every 
b G NA, IPA, c {b) is solved by a group relaxation G°{b) where a is a maximal face of the 
regular triangulation A c . 



Theorem 6.2. The following conditions are equivalent: 

(i) IPa,c ^ a Gomory family. 

(ii) The associated sets of IPa, c are precisely the maximal faces of A c . 
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(Hi) (*, t) is a standard pair of O c if and only if r is a maximal face of A c . 
(iv) All standard polytopes of IPa,c ore simplices. 

Proof. By Definition |6.1| , IPa,c is a Gomory family if and only if for all b £ NA, IPa. c {o) can 
be solved by one of its Gomory relaxations. By Theorem |3.11| , this is equivalent to saying 
that every u £ O c lies in some S(*, a) where a is a maximal face of A c and (*, a) a standard 
pair of O c . Definition [O] then implies that all associated sets of IPa, c are maximal faces 
of A c . By Theorem every maximal face of A c is an associated set of IPa,c and hence 
The equivalence of statements (ii), (iii) and (iv) follow from Theorem |3.10|. □ 



A c , then we 



If c is a generic cost vector such that for a triangulation A of cone(A), A 
say that A supports the order ideal O c and the family of integer programs IPa,c- No regular 
triangulation of the matrix A in Example |4.2| supports a Gomory family. Here is a matrix 
with a Gomory family. 

Example 6.3. Consider the 3x6 matrix 

" 1 1 1 1 1 
A= 1112 2 
1 2 3 4 

In this case, cone(A) has 14 distinct regular triangulations and 48 distinct sets O c as c varies 
among all generic cost vectors. Ten of these triangulations support Gomory families; one for 
each triangulation. For instance, if c = (0, 0, 1, 1, 0, 3), then 

A c = {a 1 = {1,2,5}, a-2 = {1,4, 5}, a 3 = {2,5,6}, a 4 = {4, 5, 6}} 

and IPa, c is a Gomory family since the standard pairs of O c are: 

(0,cxi), (e 3 ,ai), (e 4 ,ai), (0,a 2 ), (0,a 3 ), and(0,a 4 ). 

□ 

Algebraically, IPa.c is a Gomory family if and only if the initial ideal m c (/^) has no 
embedded primes and hence Theorem is a characterization of toric initial ideals without 
embedded primes. A sufficient condition for an ideal in k[xi, . . . ,x n ] to not have embedded 



primes is that it is Cohen- Macaulay JT(|. In general, Cohen-Macaulayness is not necessary 
for an ideal to be free of embedded primes. However, empirical evidence seemed to suggest 
for a while that for toric initial ideals, Cohen-Macaulayness might be equivalent to being free 
of embedded primes. A counterexample to this was found recently by Laura Matusevich. 
The algebraic approach to integer programming allows one to compute all down sets O c of 



a fixed matrix A as c varies among the set of generic cost vectors. See [24], [36] and [37] for 



details. The software package TiGERS is custom-tailored for this purpose. 

We now compare the notion of a Gomory family to the classical notion of total dual 
integrality [[33], §22]. It will be convenient to assume that /LA = Z d for these results. 

Definition 6.4. The system yA < c is totally dual integral (TDI) if LPa, c {o) has an 
integral optimal solution for each b £ cone(A) (1 Z d . 



Definition 6.5. 

face a £ A r . 



The regular triangulation A c is unimodular if ZA a = Z d for every maximal 
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Example 6.6. The regular triangulation in Example 
Example |2.2| (ii) and (iii) are not. 



21 (i 



is unimodular while those in 

□ 



Lemma 6.7. The system yA < c is TDI if and only if the regular triangulation A c is 
unimodular. 



Proof. The regular triangulation A c is the normal fan of P c by Proposition |2.4| , and it is 
unimodular if and only if 1A a = Z, d for every maximal face a G A c . This is equivalent to 
every b G cone(A a ) H Z d lying in NA a for every maximal face a of A c . By Lemma [2.3| , this 
happens if and only if LP^ c (b) has an integral optimum for all b G cone(A) n If 1 . □ 

Corollary 8.4 in shows that A c is unimodular if and only if the monomial ideal in c (Ij\) 
is generated by square-free monomials. Hence, by computing in c (/^), one can determine 
whether yA < c is TDI. Such computations can be carried out on computer algebra systems 
like CoCoA U or MACAULAY 2 |18| for moderately sized examples. See for algorithms. 



Standard pair decompositions of monomial ideals can be computed with MACAULAY 2 [pO 
Theorem 6.8. If yA < c is TDI then IPa, c is a Gomory family. 

Proof. By Theorem |4.3| , (0, a) is a standard pair of O c for every maximal face a of A c . 
Lemma |6.7| implies that cone{A a ) is unimodular (i.e., ZA a = Z d ), and therefore NA a = 
cone(A a ) fl 7L d for every maximal face a of A c . Hence the semigroups NA a arising from the 
standard pairs (0, a) as o varies over the maximal faces of A c cover NA. Therefore the only 
standard pairs of O c are (0, a) as a varies over the maximal faces of A c . The result then 
follows from Theorem 16. 2j. □ 



When yA < c is TDI, the multiplicity of a maximal face of A c in O c is one (from Theo- 
rem f4.5| ). By Theorem |6.8| , no lower dimensional face of A c is associated to IPa,c- While this 
is sufficient for IP^db) to be a Gomory family, it is far from necessary. TDI-ness guarantees 
local integrality in the sense that LPA, c {b) has an integral optimum for every integral b in 
cone(A). In contrast, if IPa,c is a Gomory family, the linear optima of the programs in LPa, c 
may not be integral. 

If A is unimodular (i.e., 7jA a = 7L d for every nonsingular maximal submatrix A G of A), 
then the feasible regions of the linear programs in LPa, c have integral vertices for each b G 
cone(A) nZ d , and yA < c is TDI for all c. Hence if A is unimodular, then IPa, c is a Gomory 
family for all generic cost vectors c. However, just as integrality of the optimal solutions of 
programs in LPa )C is not necessary for IPa, c to be a Gomory family, unimodularity of A is 
not necessary for IPa,c to be a Gomory family for all c. 

Example 6.9. Consider the seven by twelve integer matrix 



.4 



■ 1 
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of rank seven. The maximal minors of A have absolute values zero, one and two and hence 
A is not unimodular. This matrix has 376 distinct regular triangulations supporting 418 
distinct order ideals O c (computed using TiGERS). In each case, the standard pairs of O c 
are indexed by just the maximal simplices of the regular triangulation A c that supports it. 
Hence IPa,c is a Gomory family for all generic c. □ 



The above discussion shows that IPa, c being a Gomory family is more general than yA < c 
being TDI. Similarly, IPa, c being a Gomory family for all generic c is more general than A 
being a unimodular matrix. 

7. Gomory Families and Hilbert Bases 

As we just saw, unimodular matrices or more generally, unimodular regular triangulations 
lead to Gomory families. A common property of unimodular matrices and matrices A such 
that cone(A) has a unimodular triangulation is that the columns of A form a Hilbert basis 
for cone(A), i.e., NA = cone(A) R 7h d (assuming TLA = Z d ). 

Definition 7.1. Adxn integer matrix A is normalif the semigroup NA equals cone(A)T\'L d . 

The reason for this (highly over used) terminology here is that if the columns of A form 
a Hilbert basis, then the zero set of the toric ideal 1a (called a toric variety) is a normal 
variety. See 0, Chapter 14] for more details. We first note that if A is not normal, then 
IPa,c need not be a Gomory family for any cost vector c. 



Example 7.2. The matrix A 



is not normal since (1,2)* which lies in 



1111 
13 4 

cone(A) R Z 2 cannot be written as a non-negative integer combination of the columns of A. 
This matrix gives rise to 10 distinct order ideals O c supported on its four regular triangu- 
lations {{1, 4}}, {{1, 2}, {2, 4}}, {{1, 3}, {3, 4}} and {{1, 2}, {2, 3}, {3, 4}}. Each O c has at 
least one standard pair that is indexed by a lower dimensional face of A c . 



The matrix in Example [O is also not normal and has no Gomory families. While we do 



not know whether normality of A is sufficient for the existence of a generic cost vector c such 
that IPa, c is a Gomory family, we will now show that under certain additional conditions, 
normal matrices do give rise to Gomory families. 

Definition 7.3. Adxn integer matrix A is A-normal if cone(A) has a triangulation A 
such that for every maximal face a G A, the columns of A in cone(A a ) form a Hilbert basis. 



Remark 7.4. If A is A-normal for some triangulation A, then it is normal. To see this note 
that every lattice point in cone(A) lies in cone(A a ) for some maximal face a G A. Since A 
is A-normal, this lattice point also lies in the semigroup generated by the columns of A in 
cone(A a ) and hence in NA. 

Observe that A is A-normal with respect to all the unimodular triangulations of cone(A). 
Hence triangulations A with respect to which A is A-normal generalize unimodular trian- 
gulations of cone(A). 
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Examples [7.5| and |7.6| show that the set of matrices where cone(A) has a unimodular 
triangulation is a proper subset of the set of A-normal matrices which in turn is a proper 
subset of the set of normal matrices. 

Example 7.5. Examples of normal matrices with no unimodular triangulations can be 

If cone(A) is simplicial for such a matrix, A will be A-normal with 



found in flq] and [11 



respect to its coarsest (regular) triangulation A consisting of the single maximal face with 



support cone(A). For instance, consider the following example taken from [TT 



.4 



1 







Here cone(A) has 77 regular triangulations and no unimodular triangulations. Since cone(A) 
is simplicial, A is A-normal with respect to its coarsest regular triangulation {{1,2,3,8}}. 

Example 7.6. There are normal matrices A that are not A-normal with respect to any 
triangulation of cone(A). To see such an example, consider the following modification of the 
matrix in Example [75| that appears in 0, Example 13.17] : 



A 



10 11111 
1 1 1 2 2 2 
1 1 2 2 3 3 
000012345 
111111111 



This matrix is again normal and each of its nine columns generate an extreme ray of cone(A). 
Hence the only way for this matrix to be A-normal for some A would be if A is a unimodular 
triangulation of cone(A). However, there are no unimodular triangulations in this example. 



Theorem 7.7. If A is A-normal for some regular triangulation A then there exists a generic 
cost vector c £ Z n such that A = A c and IPa, c is a Gomory family. 

Proof. Without loss of generality we can assume that the columns of A in cone(A a ) form a 
minimal Hilbert basis for every maximal face a of A. If there were a redundant element, the 
smaller matrix obtained by removing this column from A would still be A-normal. 

For a maximal face a £ A, let oi n C {1, . . . , n} be the set of indices of all columns of 
A lying in cone{A a ) that are different from the columns of A a . Suppose a^,... , a ik are 
the columns of A that generate the one dimensional faces of A, and d £ W 1 a cost vector 
such that A = A c >. We modify d to obtain a new cost vector c £ MJ 1 such that A = A c 



as follows. For j = 1, 



then dj 



E 



£ R d+1 lies in C a := 
of C = cone({a\, c'i) : % — 1, 



(dj, Cj t 



.-o-AiOj, < Aj < 1 and we define Cj 



cone{{a\, : 
. . ,n). If y £ 



i £ a 

1®d 



If j £ a i n for some maximal face a £ A, 
Y.iea^i- Hence, for all j £ a in , 
cone((al, d \) : i £ a) which was a facet 



is a vector as in Definition 2.1 showing that a 



is a maximal face of A c / then y ■ ai = Ci for all i £ a U o- in and y ■ aj < Cj otherwise. Since 
cone(A a ) = cone(A aUa . m ) , we conclude that cone(A a ) is a maximal face of A c . 
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If b G NA lies in cone(A a ) for a maximal face a G A c , then IP^ )C (b) has at least one feasible 
solution u with support in cr U crj n since A is A-normal. Further, (&*, c • w) = ((Am)*, c • u) 
lies in C CT and all feasible solutions of IPA, c {b) with support in cr U Oi n have the same cost 
value by construction. Suppose v G N™ is any feasible solution of JPa jC (6) with support not 
in a U a in . Then c ■ w < c ■ t> since (a*, q) G C CT if and only if i G a U a iri and is a lower 
facet of C. Hence the optimal solutions of 7P^ iC (6) are precisely those feasible solutions with 
support in cr U a in . The vector b can be expressed as b = b' + J2i €a z i a i where Zi G N are 
unique and b' G {^ i6(T ^* ai : < A.j < 1} fl Z d is also unique. The vector b' = J2jea in r j a j 
where Tj G N. Setting U{ = Zi for all i G cr, «j = fj for all j G cr in and = otherwise, we 
obtain all feasible solutions u of IPa,c(P) with support in a U cr iri . 

If there is more than one such feasible solution, then c is not generic. In this case, we 
can perturb c to a generic cost vector c" = c + eu by choosing 1 ^> e > 0, Uj <^ 
whenever j = i\, . . . ,ik and Uj = otherwise. Suppose Ui, . . . ,u t are the optimal solutions 
of the integer programs IPA, c "(b') where b' G {^ igo . \&i '■ < A.; < 1} PI Z d . (Note that 
t = |{X^s<7 ^i a i '■ < A.j < 1} n Z d | is the index of ZA a in ZA) The support of each such 
Ui is contained in <Tj„. For any b G cone(A a ) fl Z d , the optimal solution of IPA, c "{b) is hence 
u = Ui + z for some i G {1, . . . , t} and z G N n with support in cr. This shows that NA is 
covered by the affine semigroups 4>A{S{ui, a)) where cr is a maximal face of A and Ui as above 
for each a. By construction, the corresponding admissible pairs (ui, a) are all standard for 
O c ii. Since all data is integral, c" G Q n and hence can be scaled to lie in Z n . Renaming c" 
as c, we conclude that IPa,c is a Gomory family. □ 

Corollary 7.8. Let A be a normal matrix such that cone(A) is simplicial, and let A be the 
coarsest triangulation whose single maximal face has support cone(A). Then there exists a 
cost vector c G Z n such that A = A c and IPa, c is a Gomory family. 

Example 7.9. Consider the normal matrix in Example |6.3| . Here cone(A) is generated by 
the first, second and sixth columns of A and hence A is A-normal with respect to the regular 
triangulation {{1, 2, 6}}. There are 13 distinct sets O c supported on A. Among the 13 
corresponding families of integer programs, only one is a Gomory family. A representative 
cost vector for this IPa, c is c = (0, 0, 4, 4, 1, 0). The standard pair decomposition of O c is the 
one constructed in Theorem |7.7| . The affine semigroups S(-,a) from this decomposition are: 

S(0, cr), 5(e 3 ,cr), S(e 4 ,cr), and S(e 5 , a). 

Note that A is not A-normal with respect to the regular triangulation supporting the Go- 
mory family IPa, c m Example |0|. The columns of A in cone(A ai ) are the columns of A ai 
and A3- The vector (1, 2, 2) is in the minimal Hilbert basis of cone(A ai ) but is not a column 
of A This example shows that a regular triangulation A of cone(A) can support a Gomory 
family even if A is not A-normal. The Gomory families in Theorem |?T7] have a very special 
standard pair decomposition. □ 

Problem 7.10. If A G Z dxn is a normal matrix, does there exist a generic cost vector c G Z n 
such that IPa, c is a Gomory family? 



32 



REKHA R. THOMAS 



While we do not know the answer to this question, we will now show that stronger results 
are possible for small values of d. 

Theorem 7.11. If A G Z dxn is a normal matrix and d < 3, then there exists a generic cost 
vector cGZ" such that IPa, c is a Gomory family. 

Proof. It is known that if d < 3 then cone(A) has a regular unimodular triangulation A c 
|34|1 . The result then follows from Corollary |6.8| . □ 



Before we proceed, we rephrase Problem [7.101 in terms of covering properties of cone(A) 



and NA along the lines of @, ||, 0, |TT| and [j34|- To obtain the same set up as in 



these papers we assume in this section that A is normal and the columns of A form the 
unique minimal Hilbert basis of cone(A). Using the terminology in J7|, the free Hilbert cover 
problem asks whether there exists a covering of NA by semigroups NA T where the columns 
of A T are linearly independent. The unimodular Hilbert cover problem asks whether cone(A) 
can be covered by full dimensional unimodular subcones cone(A T ) (i.e., "ZA T = Z d ), while 
the stronger unimodular Hilbert partition problem asks whether cone(A) has a unimodular 
triangulation. (Note that if cone(A) has a unimodular Hilbert cover or partition using 
subcones cone(A T ), then NA is covered by the semigroups NA r .) All these problems have 
positive answers if d < 3 since cone(A) admits a unimodular Hilbert partition in this case [||, 
||34j| . Normal matrices (with d — 4) such that cone(A) has no unimodular Hilbert partition 
can be found in M and |11|| . Examples (with d — 6) that admit no free Hilbert cover and 
hence no unimodular Hilbert cover can be found in M and H. 

When yA < c is TDI, the standard pair decomposition of NA induced by c gives a 
unimodular Hilbert partition of cone(A) by Theorem |6.7| . An important difference between 



Problem |7.10| and the Hilbert cover problems is that affine semigroups cannot be used in 
Hilbert covers. Moreover, affine semigroups that are allowed in standard pair decompositions 
come from integer programming. If there are no restrictions on the affine semigroups that 
can be used in a cover, NA can always be covered by full dimensional affine semigroups: 
for any triangulation A of cone(A) with maximal subcones cone{A,j), the affine semigroups 
b + NA a cover NA as b varies in {Xliecr^« a « : < Aj < 1} fl Z d and a varies among the 
maximal faces of the triangulation. A partition of NA derived from this idea can be found 
in 



35], Theorem 5.2]. We recall the notion of supernormality introduced in (19 



Definition 7.12. A matrix A G Z dxn is supernormal if for every submatrix A' of A, the 
columns of A that lie in cone(A') form a Hilbert basis for cone(A'). 



Proposition 7.13. For A G Z dxn , the following are equivalent: 

(i) A is supernormal, 

(ii) A is A-normal for every regular triangulation A of cone(A), 

(iii) Every triangulation of cone(A) in which all columns of A generate one dimensional 
faces is unimodular. 

Proof. The equivalence of (i) and (iii) was established in |19|, Proposition 3.1]. Definition 
|7.12| shows that (i) =^ (ii). Hence we just need to show that (ii) =>- (i). Suppose that A is 
A-normal for every regular triangulation of cone(A). In order to show that A is supernormal 
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we only need to check submatrices A' where the dimension of cone(A') is d. Choose a cost 
vector c with q ^> if the zth column of A does not generate an extreme ray of cone(A'), 
and Ci = otherwise. This gives a polyhedral subdivision of cone(A) in which cone(A') is a 
maximal face. There are standard procedures that will refine this subdivision to a regular 
triangulation A of cone(A). Let T be the set of maximal faces a of A such that cone(A a ) 
lies in cone(A'). Since A is A-normal, the columns of A that lie in cone(A a ) form a Hilbert 
basis for cone(A a ) for each a G T. However, since their union is the set of columns of A that 
lie in cone(A'), this union forms a Hilbert basis for cone(A'). □ 

It is easy to catalog all A-normal and supernormal matrices, of the type considered in this 
paper, for small values of d. We say that the matrix A is graded if its columns span an afline 
hyperplane in M. d . If d = 1, cone(A) has n triangulations {{i}} each of which has the unique 
maximal subcone cone(cii) whose support is cone(A). If we assume that a\ < a<i < • • • < a n , 
then A is normal if and only if either a% = 1, or a n = — 1. Also, A is normal if and only if 
it is supernormal. If d = 2 and the columns of A are ordered counterclockwise around the 
origin, then A is normal if and only if det(a,i, a, + i) = 1 for all % = 1, . . . , n — 1. Such an 
A is supernormal since it is A-normal for every triangulation A — the Hilbert basis of a 
maximal subcone of A is precisely the set of columns of A in that subcone. If d = 3 then 
as mentioned before, cone(A) has a unimodular triangulation with respect to which A is 
A-normal. However, not every such A needs to be supernormal: the matrix in Example |6.3| 
is not A-normal for the A supporting the Gomory family in that example. If d = 3 and A 
is graded, then without loss of generality we can assume that the columns of A span the 
hyperplane x\ = 1. If A is normal as well, then its columns are precisely all the lattice points 
in the convex hull of A. Conversely, every graded normal A with d = 3 arises this way - 
its columns are all the lattice points in a polygon in R 2 with integer vertices. In particular, 
every triangulation of cone(A) that uses all the columns of A is unimodular. Hence, by 
Proposition |7.13|, A is supernormal, and therefore A-normal for any triangulation of A. 



Theorem 7.14. Let A e Z dxn be a normal matrix of rank d. 

(i) If d = 1,2 or A is graded and d = 3, every regular triangulation of cone(A) supports at 
least one Gomory family. 

(ii) If d = 2 and A is graded, every regular triangulation of cone(A) supports exactly one 
Gomory family. 

(iii) If d = 3 and A is not graded, or if d = 4 and A is graded, then not all regular 
triangulations of cone(A) may support a Gomory family. In particular, A may not be 
A-normal with respect to every regular triangulation. 

Proof, (i) lid = 1, 2 or A is graded and d = 3, A is supernormal and hence by Proposition |7.13| 
and Theorem [7.7| , every regular triangulation of cone(A) supports at least one Gomory family. 

(ii) If d = 2 and A is graded, then we may assume that 



A 



111... 1 
1 2 ... n-l 



In this case, A is supernormal and hence every regular triangulation A of cone(A) supports 
a Gomory family by Theorem [7.7| . Suppose the maximal cones of A, in counter-clockwise 
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order, are C\, . . . ,C r . Assume the columns of A are labeled such that Cj = cone(aj_i, ai) 
for i = 1, . . . , r, and the columns of A in the interior of Cj are labeled in counter-clockwise 
order as bn, . . . , 6^. Hence the n columns of A from left to right are: 

&cb &n> " " ■ j ^lfen a ii &2i> ' ' ' , o r _i, 6 r i, ■ • ■ , b r k r , a r . 

Indexing the columns of A by their labels, the maximal faces of A are oi = {i — for 
% = 1, . . . , r. Let Ci be the unit vector of W 1 indexed by the true column index of dj in A 
and tij be the unit vector of W 1 indexed by the true column index of 6^ in A. Since the 
columns of A form a minimal Hilbert basis of cone(A), d is the unique solution to IPa,c{(H) 
for all c and e^- is the unique solution to IPA, c {bij) for all c. Hence the standard pairs of 
Theorem [Tl] are (0, ai) and (e^, Uj) for i = 1, . . . , r and j = 1, . . . , fcj. 

Suppose A supports a second Gomory family IPa,u- Then every standard pair of O w is 
also of the form (*, (jj) for a, G A, and r of them are (0, a"j) for z = 1, . . . , r. The remaining 
standard pairs are of the form (eij,<Jk). To see this, consider the semigroups in NA arising 
from the standard pairs of O w . The total number of standard pairs of O c and O w are the 
same. Since the columns of A all lie on x\ = 1, no two can be covered by a semigroup 
coming from the same standard pair and none of them are covered by a semigroup (0,0",). 
We show that if (e^, Ofc) is a standard pair of O w then k — i and thus O w = O c . 

If r — 1, the standard pairs of are (0, <Ti), (en, <Ti), . . . , (e^, cri) as in Theorem [T7| . 
If r > 1, consider the last cone C r = cone(a r _i,a r ). If a r _! is the second to last column of 
A, then C r is unimodular and the semigroup from (0, ay) covers C r fl Z 2 . The subcomplex 
comprised of C\, . . . , C r _i is a regular triangulation A' of cone(A') where A' is obtained 
by dropping the last column of A. Since A' is a normal graded matrix with d = 2 and A' 
has less than r maximal cones, the standard pairs supported on A' are as in Theorem |7.7| 
by induction. If a r _i is not the second to last column of A then b r k r , the second to last 
column of A is in the Hilbert basis of C r but is not a generator of C r . So O w has a standard 
pair of the form (e r k r , cr,). If <jj 7^ oy, then the lattice point 6 r fc r + a r cannot be covered by 
the semigroup from this or any other standard pair of O w . Hence <jj = a r . By a similar 
argument, the remaining standard pairs indexed by ay are (e r ^ r -i), u r ), • • • , (e r i,a r ) along 
with (0,a r ). These are precisely the standard pairs of O c indexed by a r . Again we are 
reduced to considering the subcomplex comprised of Ci, . . . ,C r _i and by induction, the 
remaining standard pairs of O w are as in Theorem [7.7| . 

(hi) The 3x6 normal matrix A of Example |3.3| has 10 distinct Gomory families supported 
on 10 out of the 14 regular triangulations of cone(A). Furthermore, the normal matrix 

"lllllll" 
10 11110 
12 2 110 
4 3 2 1 

has 11 distinct Gomory families supported on 11 out of its 19 regular triangulations. □ 
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